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On  the  Nuclear  Magnetic  Resonance  in  Metals  and  Alloys 

by 

N.  Bloembergen  and  T.  J.  Rowland 
Division  of  Applied  Science,  Harvard  University 
Cambridge,  Massachusetts 


Abstract 


The  shift  of  the  nuclear  resonance  and  the  line  width  has  been 
measured  in  a-  and  p-tin,  thallium,  lead  and  several  adloys.  While  gray 
tin  has  no  shift  relative  to  the  resonance  in  insulating  compounds, 
white  tin  has  an  anisotropic  shift  of  0.  79  per  cent  along  the  tetragonal 
axis  and  0.74  per  cent  perpendicular  to  it.  A theory  for  the  anisotropy 
and  asymmetry  of  the  observed  line  is  given. 

Thallium  has  an  anomalously  broad  resonance  line  with  a shift  of 
1.54  per  cent.  In  thalliiun  alloys  the  shift  varies  cxmtinuously  with 
composition  in  each  phase  and  has  discontinuities  when  phase  transitions 
occur.  A maximum  shift  of  2. 1 per  cent  has  been  found  in  a TB  Mg  - alloy, 
while  the  Na  T6  structure  exhibits  a negative  shift  of  >1.0  per  cent. 

Quadrupolar  effects  are  important,  when  the  nuclear  spin  l>^'y>  The 
intensity  of  the  copper  resonance  in  pure  copper  can  be  increased  by  an> 
nealing  and  decreased  by  cold  work.  In  alloys  like  a>brass  and  copper - 
silver  the  intensity  decreases  rapidly  with  increasing  zinc  or  silver  content. 
At  25  per  cent  zinc,  the  copper  resonance  becomes  imobservable.  These 
effects  are  explained  in  terms  of  quadrupole  interaction  and  can  give  in- 
formation about  the  state  of  order  or  disorder  on  an  atomic  scale  in  the 
alloy. 


1 . Introduction 

Nuclear  magnetic  resonance  experiments  have  already  yielded  some 
interesting  information  on  the  electronic  structure  and  self -diffusion  in 
metals.  The  three  main  sources  of  information  of  nuclear  resonance  about 
the  structure  of  matter  are  the  position  of  the  resonance  line,  the  width 
and  shape,  and  the  spin-lattice  relaxation  time. 

It  was  discovered  by  Knight^  that  the  resonance  frequency  in  metals  is 
higher  than  for  nuclei  of  the  same  isotope  in  an  insulating  material  in  the 
same  magnetic  field.  The  effect  is  caused  by  the  local  field  produced  by  the 
conduction  electrons  as  was  first  sliown  by  Townes.^  The  internal  field 
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produced  at  the  poaition  of  the  nuclei  is  much  larger  than  might  be  anticipated 
from  the  electronic  paramagnetic  susceptibility,  because  of  the  partial  S-wave 
character  of  the  wave  function.  Let  the  expectation  value  for  an  electron 
with  a wave  vector  k near  the  Fermi  level,  to  be  at  the  nucleus  in  the 
metal,  be  given  by 


^ met<0> 


= C, 


4«  » (0) 

^ atom;  ' 


(1) 


Here  4*0^^^  metallic  wave  function  normalised  to  unity  over  one  atomic 

volume  V , 4<  (0)  is  the  atomic  wave  function  where  the  conduction 

o atom 

electron  is  supposed  to  be  in  an  S-etate  outside  a closed  shell.  The  hyper > 
fine  interaction  in  the  free  atom  in  an  S-state^  is  given  by 


A I 


4>  (0) 

atom 


where  the  hyperfine  structure  constant  is  proportional  to 
Ic  the  metal  only  2B  H V conduction  electrons  per  atom  have  un- 

balanced  spins,  with  magnetic  moments  parallel  to  Here  is  the 

Bohr  magnetpn  and  N(Ep)  is  the  density  of  states,  the  nximber  of  states  per 
unit  volume  and  per  vinit  energy  interval,  at  the  Fermi  level.  For  a free 
electron  gas 


V^N(E)=|Ep‘ 


In  the  metal  the  nuclear  spin  and  the  electron  spin  will  always  be  decoupled 
d\ie  to  the  strong  interaction  between  conduction  electrons.  For  a nuclear 
spin  transition  Am^  ~ 1.  interaction  energy  with  the  electron  spins 

changes  by  an  amount 

YRAK  = hAv  = 6H  V N(E)  C A (2) 

Here  C is  a suitable  average  of  over  the  Fermi  surface,  \ is  the  nuclear 
gyromagnetic  ratio.  yfl.I  = is  the  maximum  value  of  the  s -component  of  the 
nuclear  magnetic  moment.  The  relative  Knight  shift  can  thus  be  written  in 
the  form 


-AH 

«T 


^ C A V N(E) 


o 


(3) 
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If  N(£)  IS  known  from  electronic  specific  heat  or  magnetic  susceptibility 
measurements,  the  value  of  C can  be  determined  experimentally  from  the  ob- 
served Knight  shift  and  compared  with  theoretical  values.  Ore  would  eXpoct 
the  values  of  C to  be  of  the  order  of  unity,  as  the  metallic  wave  function 
near  the  Fermi  level  will  usually  still  have  a large  percentage  of  S-wave 
character  and  the  renormalization  factor  should  not  change  the  order  of 
magnitude.  The  observed  values  of  are  in  agreement  with  the  expec- 

tation and  range  from  0.026  per  cent  for  Li  to  1.5d  per  cent  for  thallium, 
the  largest  value  reported  so  far  for  any  element.  This  is,  of  course,  directly 
related  with  the  large  value  of  atomic  hyperfine  splitting  for  this  heavy 
element. 


The  relaxation  time  in  metals  has  been  discussed  theoretically  by 

Heitler  and  Teller,^  Korringa^  and  Overhauser and  has  been  studied  experi- 

8 9 10 

mentally  by  Bloembergen,  Poulis  and  Rollin.  It  is  determined  by  the  off- 
diagonal  components  of  the  nuclear  spin -electron  spin  interaction.  In  the 
fundamental  process  these  spins  flip  simultaneously  AnOj  = - Am^  “ 1. 
balance  of  energy  . - y )H  is  taken  up  by  a change  in  wave  number  of 

the  conduction  electron.  Only  electrons  near  the  Fermi  level  can  take  up 
such  a small  amount  of  energy.  The  probability  for  the  process  is  therefore 
proportional  to  kT  and  the  complete  theoretical  expression  for  the  relaxation 
time  Tj  is 


1 irkT 
T-, 


V ^ N^E) 
o 


(4) 


From  eqs.  (3)  and  (4)  follows  Korringa’s  relation  between  the  line  shift  and 
the  relaxation  time 


1 


(5) 


The  relation  eq.  (5)  is  only  approximately  fulfilled.  Apother  term  should  be 
added  to  the  right-hand  side  of  eq.  (4)  describing  the  interaction  of  the 
nuclear  spin  with  the  non-spherically  symmetrical  part  of  theeelectronic  motion. 
The  complete  Hamiltonian  for  the  nuclear  spin  interaction  with  the  conduction 
electrons  is 
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^N 


(6) 


Here  r is  the  radius  vector  of  the  electron  with  the  nucleus  at  the  origin. 

The  first  term  is  the  point  interaction  discussed  before.  The  second 
term  represents  the  interaction  of  the  nuclear  spin  with  the  electron  spin, 
the  third  term  with  the  electron  orbit,  the  last  term  is  the  quadrupole  inter- 
action, not  included  in  Korringa's  paper.  To  find  the  matrix  elements  the 
terms  of  the  Hamiltonian  must  be  integrated  over  the  multi -electron  wave 
function  of  the  metal.  In  our  approximation  we  integrate  over  a one -electron 
wave  fxmetion  and  then  sum  over  all  wave  numbers  k.  The  orbital  angu- 
lar momentum  p in  the  metal  is  quenched.  There  is  no  preference  for  a left- 
or  right-handed  motion  of  the  electrons,  and  this  interaction  will  not  con- 
tribute to  the  Knight  shift  in  a first  approximation.  The  spin-spin  inter- 
action, the  second  term  on  the  right  side  of  eq.  (6),  will  make  a contribu- 
tion in  a non-cubic  crystal.  The  effect  will  depend  on  the  angle  of  the 
crystalline  axes  and  the  external  field,  and  the  average  for  a polycrystal- 
line specimen  will  be  zero.  This  anisotropy  of  the  Knight  shift  will  be 
discussed  in  detail  in  section  5 below.  It  leads  to  a spurious  broadening 
and  asymmetry  of  the  line  in  polycrystalline  material,  which  is  proportional 
to  the  external  magnetic  field.  If  the  crystal  is  non-cubic  and  if  the 
nuclear  spin  1^1,  the  last  term  does  not  vanish.  In  polycrystalline  ma- 
terial the  quadrupole  interaction  will  give  rise  to  a broadening  of  the  line 
which  will  be  discussed  more  fully  in  section  Z.  Clearly  the  complete  quadru- 
pole interaction  also  includes  the  effect  of  nuclear  charges  and  non-con- 
duction electrons.  All  these  contributions  to  the  inhomogeneous  electric 
field  can  be  lumped  together  in  discussing  the  line  position  and  broadening. 


The  contribution  of  the  Kairviltonian  eq.  (6)  to  the  relaxation  time  has 
a different  character.  For  the  Bloch-type  ono-electron  wave  functions 
and  non-diagonal  elements  of  the  must  be  squared, 

averaged  over  the  Fermi  distribution  and  multiplied  by  the  density  of  final 
states  to  obtain  the  total  transition  probability.  The  contribution  of  the 
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last  three  terms  is  now  nonvanishing  even  in  cubic  crystals,  since  tl'.e 
individvial  wave  function  does  not  have  this  symmetry.  This  result  has 
to  be  added  to  the  right-hand  side  of  eq.  (4),  and  the  relaxation  time 
should  be  shorter  than  that  calculated  from  eq.  (5).  This  has  been  fovind 
generally  to  be  the  case.  Since  only  electrons  near  the  tail  of  the  Terxiii 
distribution  will  be  able  tc  take  up  the  small  amount  of  energy  of  the  pro- 
cesses involved,  the  transition  probability  will  remain  proportional  to  the 
absolute  temperature  and  independent  of  the  external  field,  even  if  the  com- 
plete Hamiltonian  is  used.  The  Korringa  relation  should  be  more  closely 
obeyed,  the  more  dominant  the  S-wave  character  of  the  one-electron  wave 
functions  near  the  Fermi  level. 


The  inverse  proportionality  between  T^  and  T has  been  established 
rather  satisfactorily  in  a number  of  cases,  although  a notable  exception  is 
presented  by  Cesium  at  high  temperatures.  There  is  no  confirmation  that 
T^  is  independent  of  H^.  In  fact,  a few  results  by  Poulis  indicate  otherwise. 
More  data  on  the  H^-dependence  of  T^  are  desirable.  Deviations  from  the 
theoretical  behavior  eqs.  (4)  and  (5)  could  be  attributed  to  correlation  and 
exchange  effects,  omitted  in  the  one-electron  approximation.  It  is  not  clear 
how  the  effective  density  of  "one -electron  states"  could  depend  very  strongly 
on  H or  T,  even  if  the  interaction  between  electrons  is  taken  into  account. 


The  experimentally  observed  shifts  AH  are  indeed  proportional  to  H 
and  seem  to  be  independent  of  T,  except  for  a snnall  and  rather  trivial  volume 
effect. 

The  line -width  is  determined  by  the  nuclear  dipolar  interaction  and 
by  the  finite  life-time  T^  caused  by  the  relaxation  mechanMm.  The  second 
moment  of  the  absorption  line  due  to  dipolar  interaction  in  a polycrystal- 
line powder  is  given  by 


1 

T 


(AH^) 


Av 


-6 


of 

(7) 


The  summation  j extends  over  identical  nuclei  in  the  crystal,  the  summation  f 
includes  all  different  nuclei.  Gutowski^^  has  shown  that  in  metals  at  higher 
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temperatures  the  contribution  of  the  relaxation  mechanism  to  the  width 
amounting  to  is  important.  In  that  case  the  width  will  increase 

proportionally  with  Ihe  temperature.  The  relaxation  mechanism  will  nearly 
always  be  the  width-determining  factor  above  the  melting  point.  In  the  liquid 
state  the  dipolar  interaction  of  the  nuclei  has  to  be  averaged  over  the  Brownian 
motion  and  is  reduced  to  a very  small  value.  A decrease  in  the  dipolar  width 
has  been  observed  with  increasing  temperature  in  the  solid  metals  lithium  and 
sodium,  when  the  self -diffusion  becomes  sufficiently  rapid.  This 

phenomenon  is  entirely  analogous  to  the  motional  narrowing  occurring  in 
liquids  and  could  be  expected  to  occur  in  other  metals  just  below  the  melting 
point.  It  gives  information  about  the  self -diffusion  as  a function  of  tempera- 
ture when  the  elementary  jumping  process  of  the  diffusion  has  a frequency 
larger  than  T*  ' . 

It  seems  profitable  to  extend  the  observations  of  the  shift,  width  and 
relaxation  time  to  alloys,  in  order  to  gather  information  about  the  density  of 
states  and  electronic  structure  in  these  systems.  In  the  next  section. we 
shall  discuss  the  importance  of  quadrupole  interaction  in  alloy  systems  and 
the  difficulties  presented  by  it. 

The  equipment  consisted  of  a Pound-Knight- Watkins  type  of  radiofre- 
14  15 

quency  spectrometer  ’ and  a permanent  magnet  with  a field  near  6000 
gauss  in  a 1-^inch  gap,  the  diameter  of  the  pole  faces  being  4 inches.  The 
field  was  calibrated  by  the  proton  or  La  resonance  in  an.  aqueous  solution 
of  Li  CB . The  inhomogeneity  in  the  volume  occupied  by  the  samples  was  less 
than  0.4  gauss ; or  negligible  compared  to  the  width  of  most  resonances  de- 
scribed in  this  paper.  The  field  in  the  gap  was  found  to  vary  with  temperature, 
about  -1.3  gauss  per^C,  probably  due  to  the  expansion  of  the  magnet  yoke.  The 
temperature  of  the  magnet  was  always  recorded  and  a correction  to. the  mag- 
netic field  was  applied.  A few  measurements  were  made  at  fields  between 
2000  and  5500  gauss  with  an  electromagnet.  This  magnetic  field  was  less 
homogeneous  and  drifted  with  time.  Further  measurements  at  other  field 
strengths  will  be  carried  out  with  an  improved  electromagnet.  The  frequency 
of  the  oscillator  was  compared  with  the  frequency  of  a standard  frequency 
meter  with  an  estimated  error  of  less  than  100  cps. 
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The  magnetic  field  was  modulated  sinusoidally  at  280  cps  with  ampli- 
tudes ranging  between  0.  2 and  20  gauss.  To  detect  the  weak  broad  lines  the 
larger  modulation  amplitudes  are  necessary.  The  frequency  of  the  oscillator 
was  slowly  swept  through  the  resonance  and  its  output  fed  into  a phase-sensi- 
tive detector  with  a time  constant  of  8 sec  (effective  band-width  about  -^cps). 
The  derivative  of  the  absorption  in  the  sample  as  a function  of  frequency  was 
recorded  on  an  Esterline -Angus  recording  instrument. 

Due  care  should  be  taken  that  the  experimentally  observed  quantities 
really  correspond  only  to  nuclear  absorption.  This  will  be  true  only  when 
the  metallic  particle  sixe  is  small  compared  to  the  skin  depth.  If  this 
condition  is  not  fulfilled,  part  of  the  observed  signal  may  be  due  to  eddy- 
current  losses  and  corrections  should  be  applied  in  order  to  evaluate  the 
correct  line  shift  and  relaxation  time.  The  metallic  samples  are  always 
used  in  a finely  divided  polycrystalline  form  to  increase  the  effective 
volume  as  the  radiofrequency  field  cannot  penetrate  further  than  a few  skin 
depths.  Although  experiments  on  single  crystals  would  be  desirable  it  would 
be  difficult  to  make  the  effective  volume  of  the  sample  sufficiently  large. 

The  samples  were  made  from  powdered  metals  suspended  in  mineral 
oil.  Oxygen-free  electrolytic  copper,  thallium,  lead,  tin  (c.p.  grade), 
mercury  (triply  distilled),  and  sine  (s.p.)  were  obtained  commercially. 
Alloys  were  made  by  melting  weighed  ratios  of  the  constituents  in  an  evac- 
uated quartz  or  glass  tube.  They  were  allowed  to  cool  slowly  and  were  held 
at  appropriate  annealing  temperatures  to  bring  about  equilibrium  in  the 
phases  under  consideration.  The  metals  were  then  filed  and  sieved  through 
a 200-mesh  sieve.  The  particle  size  was  smaller  than  the  skin  depth  in  all 
cases,  and  it  could  be  assumed  that  the  whole  weight  of  the  powdered  sample 
contributed  equally  to  the  absorption  signal.  Insulation  between  the  particles 
was  achieved  either  by  a very  thin  oxide  coating  or  by  tlie  oil  in  which  they 
were  suspended.  In  some  instances  the  powders  were  annealed  in  an  inert 
atmosphere.  In  the  oil  suspension,  particles  could  even  be  heated  above 
their  melting  point  without  agglomeration.  The  glass  tube  containing  the 
suspension  was  inserted  in  a shielded  radiofrequency  coil.  The  assembly 
could  be  inserted  in  a Dewar  flask  with  a tail  of  10-mm  inside  aim  18-mm 
outside  diameter  for  measurements  at  77°K.  In  another  assembly  a heater 


TR186 


-8- 


wire  waa  wrapped  aroxind  the  ahield  of  the  coil.  A thermal  insulating  layer 
on  the  outside  made  it  a miniature  oven.  The  temperature  could  be  raised 
to  350°C.  The  modulation  coils  were  outside  the  Dewar  or  oven.  The  coil 
was  connected  to  the  oscillator  by  a coaxial  cable  20  inches  long. 

2.  Quadrupole  Effects 

When  it  was  attempted  to  detect  the  magnetic  resonance  of  the  copper 
isotopes  in  brass  (80%  Cu  and  20%Zn),  the  signal  was  unobservable,  although 
a similar  volume  of  pure  copper  gave  a signal  to  noise  ratio  of  over  one 
hundred.  This  rather  unexpected  result  can  be  explained  by  considering 
quadrupole  interaction.  Although  the  crystallographic  structure  of  a-brass 
as  determined  by  X-rays  has  cubic  synrunetry,  the  environment  of  individual 
copper  nuclei  in  the  solid  solution  will  not  have  this  symmetry.  If  a few  of 
the  nearest  neighbors  are  Zn-atoms  the  cubic  symmetry  has  certainly  dis- 
appeared on  a microscopic  scale.  There  may  then  exist  gradients  of  the 
electric  field  which  will  vary  in  direction  and  magnitude  at  the  position  of 
the  various  copper  nuclei.  If  a copper  atom  happens  to  have  twelve  other 
copper  atoms  as  nearest  neighbors,  there  may  still  be  a deviation  from  cubic 
symmetry  arising  from  Zn  atoms  among  the  next  nearest  or  further  neigh- 
bors, adthough  in  this  case  the  magnitude  of  the  gradient  can  be  expected  to  be 
smaller.  Even  in  pure  copper  a certain  inhomogeheity  of  the  internal  electric 

field  will  exist  due  to  strains  around  dislocations  or  other  imperfections  in 

15  18 

the  lattice.  It  has  been  shown  by  Pound  and  Watkins  ’ that  these  strains 
in  single  crystals  of  the  alkali  halides  are  sufficient  to  wipe  out  the  satel- 
lite transitions,  i.e.  , the  resonance  of  bromine  and  iodine  nuclei  in  KBr 
and  K1  crystals  is  due  solely  to  transitions  between  the  levels  with  mj=-^ 
and  mj  = The  reason  is  that  the  energy  difference  between  these  two 

levels  is  unaffected  by  the  quadrupole  interaction,  as  the  first-order  approxi- 
mation quite  generally  depends  on  m^  only.  There  is  a first-order  pertur- 
bation of  the  other  energy  differences,  and  these  transitions  are  spread  out 
over  such  a large  frequency  interval  as  to  become  unobservable.  It  is 
possible  that  in  pure  and  well -annealed  metals  the  quadrupole  broadening 
is  considerably  less  and  that  the  satellite  transitions  also  contribute  to  the 
signal.  Evidence  for  this  will  be  presented  below.  It  is  true  that  the  metals 


Figure  1.  A recording  of  the  derivative  of  the  absorption  of  the  Cu 
resonance  at  roonr^  temperature  in  upper:  pure  annealed  copper;  middle: 
99%  Cu  and  1%  Ag.  [Note  the  appearance  of  weak  tails,  absent  in  pure 
copper.]  lower:  80%  Cu  - 20%  Zn.  An  extremely  weak  resonance  wUh  an 
intensity  less  than  1%  of  that  in  pure  copper  has  been  observed  at  77  K in 
this  alloy. 
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will  have  at  least  the  same  density  of  dislocations  as  the  alkali  halide 
crystals.  But  the  electron  orbits  in  the  metal  ions  are  probably  deformed 
less  under  the  influence  of  strains  than  of  the  Br~  and  I~  ions.  Further- 
more,  the  q^drupole  moments  of  the  Cu  isotopes  are  a factor  2 smaller  than 
that  of  Br.  Finally,  the  effect  of  strains  in  metals  on  the  electric  field 
gradient  at  a nucleus  a interatomic  distances  away  will  be  smaller 
because  atoms  are  displaced  rather  than  positive  or  negative  ions  carry- 
ing a net  charge. 

The  disappearance  of  the  resonance  line  in  alloys  with  a spin 
1 = -Ti  etc. , will  then  occur  in  two  steps.  At  first  the  quadrupole  inter - 
action  will  wipe  out  all  transitions  except  the  m = transition  by  a 

first-order  perturbation  process.  We  shall  write  down  the  formulas  for  the 
case  of  axial  symmetry  of  the  inhomogeneous  field.  This  will  be  correct 
if  one  salute  atom  is  in  particular  responsible  for  the  field.  For  somewhat 
higher  solute  concentrations  in  the  solid  solutions  it  is  by  no  means  correct. 
Nevertheless,  the  right  order  of  magnitude  of  the  line  displacement  and 
broadening  will  be  obtained  in  this  case.  Let  0 be  the  angle  between  the 
external  magnetic  field  and  the  axial  electric  field.  The  first-order  quad- 
rupole perturbation  gives  for  the  frequency  of  the  transition 

“ ’o  ■ (8) 

It  is  seen  that  the  frequency  of  the  transition  m = -^  is  not  changed 

in  this  approximation.  The  satellite  lines  are  broadened  since  all  angles 
0 will  occur  in  the  polycrystalline  sample.  Summing  over  all  angles  we 
find  that  each  satellite  line  has  a second  moment  in  a polycrystalline  powder 


(2m-l)^ 


1^21-1)^^ 


4 


(9) 


On  the  average  these  lines  are  not  displaced  since  ^1  - 3cos  = 0 . 

The  line  shape  is  quite  characteristic  and  has  been  described  in  other  papers. 
Often  the  broadening,  of  which  eq.  (9)  gives  the  magnitude,  will  make  the 
lines  iinobservable.  The  transition  m = 7 broadened  in 
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a second  appiroximation  according  to  the  equation 


aWj  j = h + a (1  - 9coa^0)  (1  - coe^O) 

I 


A 


21  3 

41^(21-1) 


4^2  2 
e Q q 


(10) 


(H) 


The  center  of  gravity  of  the  line  will  be  displaced  toward  lower  frequencies 
by  an  amount 


iv  — Ah 


tr/2 


-9cos^0)  ( 1 -cos^O)  sinO  dO 


-^-Ah  (12) 


while  the  second  moment  of  line  averaged  over  all  angles  is 


(6»  + ^ ll-  = 0.93  A^h"^ 


(13) 


The  experiment  in  brass  (20%  Zn)  was  carried  out  at  a frequency  = 6.3 
Mc/s.  The  line  should  have  been  at  least  50  kc  wide  in  order  to  escape 
observation.  This  enable*  us  to  calculate  a lower  limit  for  q or  grad  E in 
the  alloy.  Substituting  the  values  = - 0.157  x 10  ^^cm^  (or 

= - 0.  145  X 10  ^^cm^)  and  ^ = 50  kc  in  eq.  (9),  we  find  that 
e^Qq  ^ 0.  17  X 10^  cps  is  necessary  for  the  satellites  to  disappear.  A 
mhre  stringent  requirement  is  obtained  from  eq.  (13).  The  condition  that 
the  line  m = ^ ^ disappear  gives 


^ Q.a>  2.6  X 10^ 


cps 


or 

21  - 3 

q = 4.8x10  cm 

The  simplest  picture  one  can  form  about  the  origin  of  the  internal 
electric  field  is  that  the  Zn  atoms  have  a different  charge.  The  bivalent 
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Zn  atoms  will  produce  one  extra  positive  electronic  charge  at  their  lattice 
sites  compared  to  the  monovalent  Cu  atoms  The  electrons  around  the  Za 
atom  are  redistributed  around  the  Zn  atom  — either  by  the  formation  of  a 
bound  state  or  a distribution  of  conduction  electrons  — to  produce  a screened 
Coulomb  potential  of  the  form^^*^^ 

V = J (U) 

8 - 1 

where  a reasonable  value  for  K = 10  cm  Differentiating  this  expression 
twice , we  obtain  for  the  gradient 

8E 

eq  = -~  = e‘*^'  ( 1 + Kr  + ^ K^r^)  (15) 

One  Zn  atom  would  produce  a gradient  at  the  nearest  Cu  atom  (r  = 2.55x 
-8  22-3 

10  cm)  of  q = 6.4  x 10  cm  . This  is  much  smaller  than  the  experimen- 
tal limiting  value.  The  effect  should,  moreover,  depend  strongly  on  the 
effective  valency  of  the  alloying  agent.  Actually,  it  is  observed  that  silver 

produces  at  least  as  large  a quadrupole  broadening  as  sine.  Although 
22 

Huang  has  shown  that  considerable  polarisation  effects  may  occur  even 
if  the  solvent  and  solute  atom  have  the  same  valency,  a description  in  terms 
of  distorted  electron  orbitals  is  more  appropriate,  as  will  be  discussed 
below . 

It  is  of  interest  to  investigate  how  the  apparent  line  intensity  decreases 
as  the  concentration  of  Zn  or  Ag  increases.  In  Figure  2 the  maximum 
absorption,  corrected  for  the  effective  number  of  copper  nuclei  in  the  sample, 
is  plotted  as  a function  of  the  solvent  concentration  1-  C. 

In  pure  copper  samples  the  intensity  varies  between  0.4  and  1.0.  The 
unit  intensity  was  assigned  to  a well-annealed  pure  electrolytic  copper  sample. 
Culd-worked  copper  filings  then  had  an  intensity  between  0.4  and  0.5  with 
respect  to  this  reference  point.  The  intensity  could  be  increased  by  a factor 
of  nearly  two  when  the  edd-worked  powder  was  annealed.  The  strains  con- 
nected Mrith  dislocations  in  the  worked  material  produce  electric  field  gradi- 
ents which  wipe  out  the  transitions  m^  = 


3 

I 


1 

^7 


1 3 

and  ntij  = - by  fir  it- 


TR.ia6 


-12- 


order  quadrupole  perturbation.  JThe  intensity  ratio  of  these  satellites  to 
the  central  line  m = 3:4:3  . The  renmaining  intensity  of 

0.*^  is  thus  ascribed  to  the  m - transition. 

In  the  cold-worked  powdered  alloys  the  intensity  then  decreases 
further  by  second  order  quadrupole  perturbation  of  the  central  line.  How 
ever,  the  line  does  not  broaden  gradually  nor  does  it  show  a shift  toward 
lower  frequencies,  as  could  be  expected  from  eqs.  (12)  and  (13).  Only  a 
vestige  of  a broadened  tail  is  apparent  in  the  recording  for  a Cu^^g  alloy . 

We  assume,  therefore,  that  the  electric  field  gradient  is  very  large  and 
makes  the  line  unobservable,  when  a Zn  or  Ag  atom  is  very  close  to  a 
copper  atom,  but  is  rather  small  and  does  not  affect  the  central  line,  when 
the  Zn  or  Ag  atoms  are  further  away.  Let  the  Zn  atom  in  n neighboring 
positions  around  a copper  nucleus  make  its  resonance  unobservable,  and 
let  the  effect  of  the  Zn  atom  be  sero,  if  it  is  not  in  one  of  these  n positions. 
We  calculate  the  probability  that  a copper  atom  has  no  Zn  atom  in  its  neigh- 
boring lattice  sites.  If  the  relative  concentration  of  the  solute  is  C.-  each 
lattice  site  in  the  random  solute  solution  of  a-brass  has  a probability  1-C 
to  be  occupied  by  a Cu  atom.  The  desired  probability  is  therefore  (1-C)°. 

In  this  case  the  copper  atom  is  surrounded  only  by  copper  atoms.  The 
immediate  neighborhood  has  cubic  symmetry  and  the  copper  nucleus  will 
contribute  to  the  resonance  signal.  There  are,  of  course,  a few  other  con- 
figurations with  cubic  symmetry  (e.g.^all  neighbors  are:Zn  atoms)  but  for 
low  Zn  concentrations  these  have  a much  smaller  probability  and  need  not 
be  considered  here.  We  therefore  find  for  the  apparent  intensity 

6n27^  = n tn  ( 1 - C)  . 

In  Figure  3 we  have  plotted 7 versus  the  composition  on  a double  logarithmic 
scale.  It  is  seen  that  the  experimental  points  fall  suprisingly  well  on  a 
straight  line.  Its  slope  determines  n = 18.  If  this  number  came  out  much 
larger  ou-  theory  would  be  incorrect  as  the  underlying  assumption  that  was 
only  relatively  few  neighbors  produce  an  appreciable  gradient.  The  face- 
centered  cubic  laUice  of  copper  has  12  nearest  neighbors  at  a distance  of 
2.55  A and  six  next  nearest  neighbors  at  a distance  3.  61  The  number  of 
these  positions  agrees  well  with  the  experimental  value.  The  quadrupole 


iOO  95  90  85  80%  Cu 

6 3 

Figure  2 The  intensity  of  the  Cu  magnetic  resonance  in  annealed  and 
cold  worked  copper  and  in  alloys  with  zinc  and  silver.  By  cold  work  the 
contribution  of  the  transitions  m.  = 2 washed  out. 

Annealing  restores  the  total  intensity  of  the  nuclear  resonance. 
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interaction  in  one  of  these  18  positions  should  then  be  at  least  2.6  x iCv  cps 

as  derived  from  eq.  (13).  This  corresponds  to  an  asymmetry  in  the  wave 

9 2 2 

function  around  the  Cu  atom  equivalent  to  4 per  cent  an  atonuc  3d  4s  D 


— state  or  20  per  cent  of  the  atomic  3d^^  4p 


3/2 


state . 


23 


3/2 


In  the  case  of 


a neighboring  silver  atom  such  a distortion  is  probably  not  unreasonable.  In 

the  case  of  a neighboring  Zn  atom  the  extra  electron  may  be  shared  by  the 

surrounding  copper  atoms,  each  taking  part  with  an  unsaturated  covalent  bond. 

An  excess  of  a 20  per  cent  4p  bond  in  the  direction  of  the  Zn  atom  may  not  be 

25 

unreasonable.  Townes  and  Oailey  have  shown  that  covalent  bonding  can 
satisfactorily  describe  quadrupole  coupling  in  molecules.  We  leave  the 
question  open  of  what  description  of  distorted  orbits  will  be  most  appropriate 
for  the  metals.  Orbital  distortion,  however,  seems  to  be  the  dominant  origin 
for  quadrupole  effects  even  in  alkali  halide  crystals,  strained  by  the  presence 
of  dislocations.  Watkins  has  shown  that  a calculation  based  on  a distribution 
of  point  charges  deviating  from  cubic  symmetry  does  not  give  the  correct 
answer.  A phenomenological  piultiplication  factor ^ , larger  than  unity,  has 
to  be  introduced  to  obtain  results  in  agreement  with  the  experimental  quad- 
rupole coupling. 


24 


It  can  be  expected  that  the  shells  of  the  alkali  ions  will  be  much  less 

susceptible  to  deformation  than  the  shells  of  the  transition  elements  or 

heavier  elements  oiu  the  right-hand  side  of  the  periodic  System.  In  general, 

the  quadrupole  interaction  provides  a serious  handicap  in  the  investigation 

of  alloy  systems  with  nuclear  spin  The  line -broadening  may,  however, 

not  be  prohibitive  for  isotopes  with  small  values  of  the  quadrupole  inter- 
6 7 

action  like  La  and  Li  . In  Li -systems  the  quadrupole  effect  could  perhaps 

be  studied  in  detail  as  the  resonance  lines  may  net  disappear.  The  intensity 

of  the  nuclear  resonance  in  alloys,  if  gives  information  about  the  amount 

^ 1 

of  order  in  the  immediate  neighborhood  around  nuclei  with  spin  Tenden- 

cies of  solute  atoms  to  stay  away  from  each  other,  the  first  stages  of  precipi- 
tation and  deviations  from  perfect  solid  solution  or  perfect  order  could  be 
studied  in  this  fashion. 


In  ordered  structures  with  cubic  symmetry  the  quadrupole  broadening 
should  be  absent.  Attempts  to  detect  the  resonance  in  A-braes  near  the  50%  - 
50%  composition  have  failed,  presumably  because  the  ordering  is  not  complete. 
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AIbo  liquid  structures  should  be  accessible  to  investigation  as  *n  this  case 
the  quadrupole  interaction  has  to  be  averaged  over  the  rapidly  changing 
configurations  and  effectively  cancels,  although  there  will  remain  a con- 
tribution to  the  relaxation  process.  A small  quadrupole  contribution  to  the 

relaxation  process  is  possibly  responsible  for  the  fact  that  the  relaxation 

69  7 1 85  87 

times  for  the  pairs  of  isotopes  Ga  - Ga  and  Rb  - Rb  have  a somewhat 
smaller  ratio  than  the  inverse  ratio  of  the  squares  of  the  gyromagnetic  ratios 
for  these  isotopes.  It  so  happens  that  the  isotope  with  the  smaller  magnetic 
moment  has  the  larger  quadrupole  moment  in  both  cases.  For  gallium  we  have 


The  experimental  result  is  explained  if  the  quadrupole  interaction  contributes 

7 I 

9 per  cent  to  the  relaxation  mechanism  for  Ga  and  consequently  35  per  cent 
69 

for  Ga.  Since  the  two  rubidium  isotopes  have  different  total  spin,  the  com- 
parison is  not  as  straight-forward,  and  detailed  investigation  o transitions 
between  the  var;.ous  m-levels  would  have  to  be  made,  not  warranted  by  the 
experimental  uncertainties. 

The  use  of  isotopes  with  spin  I however,  avoids  all  complications 
of  the  quadrupole  effect.  The  fact  that  the  intensity  of  the  resonance  is  not 
drastically  reduced,  as  in  the  case  of  thallium  alloys  as  discussed  in  a later 
paragraph,  provides  perhaps  the  most  convincing  proof  of  the  assumption  of 
quadrupole  interaction  in  brass. 

3.  Nuclear  Resonance  in  Tin,  Thallium  and  Lead 

During  the  first  attempts  no  resonance  was  found  in  metallic  thallium. 

If  the  line  width  were  determined  by  nuclear  dipolar  interaction  the  second 
moment  of  the  line  eq.  (7)  woxild  be  4.  6 kc/s,  and  the  resonance  easily  ob- 
servable. The  hypothesis  was  made  that  the  spin-lattice  relaxation  time  would 
be  extremely  short  and  the  line  correspondingly  broader,  on  account  of  the  very 
large  hyperfine  splitting  in  this  element.  A broad  resonance  of  33  kc/s  between 
points  of  maximum  slope  (eqxiivalent  to  14  gauss)  was  then  found  at  liquid  air 
temperature.  Subsequent  investigations  showed  that  the  line  width  was  indepen- 


(n  (i-c) 


Figure  3 A doubly  logarithrriic  plot  of  the  intensity  of  the  central 
transition  line  m.  = '^■->-2°^  Cu°^  versus  zinc  concentration  C in 
in  a-brass.  The  intensity  is  normalized  to  unity  for  no  zinc  content. 
The  drawn  line  (1  - CP  has  a slope  corresponding  to  n = 18  with 
surprisingly  good  precision.  This  can  be  considered  as  a proof  that 
Zn  atom  shares  an  unsaturated  covalent  bond  with  the  twelve  nearest 
and  six  next  nearest  copper  atoms  only. 
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dent  of  tennperature  in  the  range  from  77°K  - 400°K.  This  seems  to  rule  out 
the  relaxation  mechanism  as  the  line-broadening  agent.  If  the  conduction 
electrons  are  responsible  for  the  relaxation,  the  line  width  should  be  propor- 
tional to  the  absolute  temperature.  For  a Raman  scattering  process  of  vibra- 

2 

tional  lattice  phonons  a T -dependence  would  be  expected. 

It  has  been  suggested  by  Gutowski^^  that  an  anisotropy  of  the  Knight  shift 
could  be  responsible  for  the  width  in  metals.  It  will  be  shown  in  the  next 
paragraph  that  the  anisotropy,  which  would  have  to  be  about  20  per  cent  in  this 
case,  produces  an  asymmetrical  line  with  a width  proportional  to  the  resonance 
frequency.  Experiments  in  external  fields  ranging  from  2000  to  6200  gauss 
revealed  no  change  in  the  width.  Furthermore,  the  line  should  become  narrow 
in  the  cubic  phase  of  thallium,  which  ic  stable  above  232°C.  No  resonance 
was  observable  above  150°C.  If  the  line  had  the  dipolar  width  of  4,6  kc/s, 
it  would  have  been  observed.  Even  more  striking  is  th<;  failure  to  detect 
the  resonance  in  the  liquid  phase,  the  melting  point  being  at  302°C. 

Although  the  width  makes  the  determination  of  the  line  position  some- 
what inaccurate,  the  relative  Knight  shift  is  so  large  that  it  can  be  n4..asured 
with  fair  precision.  Figure  4 shows  a record  of  line  derivatives  of  the 
TC^^^and  resonance  in  fields  of  about  3000  gauss  and  5600  gauss 

respectively.  These  resonances  were  compared  with  the  narrow  lines  in 
a saturated  aqueous  solution  of  thallium  acetate.  To  what  extent  these  lines 
can  really  serve  as  standards  is  open  to  criticism  in  view  of  the  concentration 
dependent  chemical  shift  occurring  in  solutions  of  thallium  salts.  For 

a Knight  shift  flr/v  = 1-54%  + .05%  has  been  found,  the  largest  so 

° 203 

far  reported  for  an  element.  The  behavior  of  the  Tt  resonance  is  baffling. 

In  going  from  2000  to  6000  oersted,  the  Knight  shift  changes  from  1 per  cent 

to  1 . 5 per  cent,  the  width  changes  from  17  kc/s  to  33  kc/s  and  the  intensity 

. 205 

decreases  by  a factor  three.  No  such  anomalies  are  observed  for  the  Tc 
resonance . 

The  two  isotopes  have  magnetic  moments  differing  only  by  1 per  cent 
and  it  is  not  clear  how  they  could  behave  so  differently,  if  our  interpretation 
of  the  recordings  is  at  all  correct.  The  preliminary  experimental  data  are 
mentioned  here  only  for  the  sake  of  completeness.  Some  numerical  values 
are  compiled  in  Table  I. 
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For  in  metallic  lead  a Knight  shift  of  1.24%  ^ 0.04%  has  been 

found.  A preliminary  value  was  previously  reported  by  Khight. 

The  Knight  shifts  for  the  resonances  of  Sn^^^  and  Sn^^^  in  white  tin 
are  0.74  per  cent  and  0.72  per  cent,  which  can  be  considered  equal  within 
the  experimental  error.  No  values  for  tin  and  thallium  had  been  reported 
previously,  a fact  undoubtedly  due  to  the  unexpected  broadness  of  the  lines. 


117  119 

The  resonances  of  Sn  and  Sn  in  grey  tin  reveal  no  Knight  shift. 

It  has  been  shown  that  the  diamond  structure  of  grey  tin  is  semiconducting. 

The  conduction  electrons  contribute  little  or  nothing  to  the  magnetic  suscepti- 
bility in  this  • -phase, X _ = - 0.  025  x 10  ^ and  X.  e = + 0. 025  x 10  ^ 
cgs  units/gram.  The  contribution  to  the  Knight  shift  is  correspondingly 
negligible  in  a-Sn.  The  nximber  of  electrons  and  holes  with  unbalanced  spin 
is  very  small,  and  furthermore  the  wave  function  has  predominantly  p-character 
near  the  top  of  the  valance  band  an^  the  bottom  of  the  conduction  band  in  the 
diamond  structure.  The  intensity  of  the  lines  for  the  two  isotopes  in  both 
modifications  is  about  equal,  in  agreement  with  their  abundance  ratio  of  1.08. 
The  lines  in  a-Sn  have  a normal  shape,  but  the  width  is  about  six  times  as 
large  as  the  dipolar  broadening.  This  can  clearly  not  be  blamed  on  the 
conduction  electrons  in  this  semiconducting  material. 


In  the  p- modification,  however,  both  resonances  have  a marked  asym- 
metry. This  may  be  understood  as  the  anisotropy  of  the  Knight  shift  in  the 
markedly  tetragonal  white  tin  structure.  A striding  confirmation  for  the 
correctness  of  this  hypothesis  would  b«  to  check  the  proportionality  of  this 
asymmetric  broadening  with  the  external  field.  The  integrated  line,  assuming 
that  the  recording  gives  the  derivative,  is  sketched  in  Figtire  6.  This  line 
will  be  interpreted  in  the  next  paragraph.  The  result  is  that  the  figures 
quoted  for  the  Knight  shift  above  are  valid,  if  the  field  is  perpendicular  to 
the  tetragonal  axis.  For  fields  parallel  to  the  axis  the  Knight  shifts  are 
0.79  per  cent  and  0.77  per  cent  or  7 per  cent  higher. 

4.  Anisotropy  of  the  Knight  Shift 

The  first  term  in  the  Hamiltonian  eq.  (6)  does  not  depend  on  the  di- 
rection of  the  external  field.  The  second  term  produces  a shift  of  the  nuclear 
resonance  which  depends  on  the  orientation  of  the  external  field  with  respect 
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Table  I 

External 

Field 


amole  Isotope 

Strength 
( Gau  s s) 

Resonance 
Frea.  Mc/s 

nm 

Line  yidth 
obs  . 

calc . 

1 

p-Sn 

6102 

9.761 

0.79  11 

1 

6102 

9. 323 

0.77  11 

Sn“’ 

6102 

9.756 

C.74  l_ 

3-5 

.71 

S„“’ 

6102 

9. 319 

0.72 

3-5 

a-Sn 

6100.  2 

9.6786 

0.00 

5 

.75 

Sn"'' 

6100.  2 

9.245 

0.00 

5 

Tt 

^205 

1993 

4.968 

1.43^.  10 

30^3 

4.6 

3007 

1.50 

34t2 

3665 

9.140 

1.54 

33 

4890 

1.61 

33 

5563 

1.55 

33 

6108 

1.53 

33 

6495 

1.70 

33 

Tf20  3 

1993 

4.899 

.99 

17 

3,9 

3050 

7.500 

1.04 

22 

3720 

9.140 

1.0^.  1 

28 

4960 

12.200 

1.07 

35^5 

5563 

13.686 

1.11 

31 

6490 

16.007 

1.33 

33 

Pb 

5553.4 

5.0026 
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to  the  crystallographic  axes.  It  vanishes  only  for  cubic  symmetry.  The 
third  term  in  eq.  (6)  can  be  neglected  as  the  orbital  momentum  is  frozen 
in.  In  the  first-  order  F>ertubation  calculation  we  are  only  interested  in 
elements  diagonal  in  the  nuclear  and  electron  spin.  If  we  donote  the  angle 
between  H and  the  radius  vector  between  electron  and  nucleus  by  a we 
retain  only  the  following. part  of  the  HamiUnnian 


»^anis  ^ ^ 

The  - or  -f  sign  depends  on  whether  the  electron  moment  is  parallel  or  anti- 
parallel  to  respectively.  To  calculate  the  shift  we  determine  the  energy 
difference  due  to  eq.  (16)  between  two  nuclear  spin  orientations  for  which 
= 1.  We  have  to  integrate  eq.  (16)  over  the  electron  wave  function 
and  then  sum  over  all  states  k,  taking  the  contributions  of  electrons  of  op- 
posite spin  with  opposite  sign.  Per  unit  volume,  only  H^N(Ep)  electrons 
near  the  Fermi  level  with  unbalanced  spin  will  therefore  contribute.  If 
is  the  atomic  volume  we  find  for  the  energy  difference 


AW  = ypjftpZ  V^N(Ep)<  J ij,*(3  cos^a^^-l)rj^'\ 


Average  (17) 
over  Fermi 
surface 


If  we  introduce  an  average  wave  function  4^,  so  that  i)'  4^  represents  the  average 
electron  density  in  space  of  the  conduction  electrons  near  the  Fermi  level,  we 
find  from  eq.  (17)  for  the  relative  Knight  shift 


anis 


H 


(3cos^o-l).‘  ^4<dxdydz 


(18) 


The  field  makes  the  polar  angles  0 and  4 with  respect  to  the  x,y,z 
coordinate  system.  The  radius  vector  j;  has  polar  angles  (§)and  6.  We  can 
now  express  (1-3  cos^a)  by  means  of  a well-known  addition  theorem  in  terms 
of  tesseral  harmonics  in  the  angles  0 and  4 and  @ and  i respectively. 

2 

■^(3  cos^o- 1)  = ^ (-l)*"p2”'(cos@)P2*"'(cosO)e'"'^*‘  (19) 

m=-2 
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We  ahall  first  discuss  the  case  of  axial  symmetry.  It  is  probably  a good 
approximation  to  assume  that  4*  consists  of  a mixture  of  p-type  wave  functions 
of  the  form  4»^,  ^*^+1  1^"  '*^ave  functions 

are  real  and  represent  the  quenching  of  the  p-orbits  by  the  crystalline  field. 
In  the  case  of  axial  symmetry  the  electron  density  can  then  be  written  as 


H = g(r) 


A<x^  + y^)  + C 


g(r)  [A  + (C  - A)cos^@  ] 


(20) 


Here  g(r)  is  a radial  function,  whose  form  is  unimportant  for  our  purposes. 
Clearly  eq.  (20)  could  be  extended  to  higher-order  terms  with  axial  symmetry, 
and  the  following  integrations  could  easily  be  carried  through,  if  d - or 
higher  orbital  states  were  included.  SubstUuting  eqs.  (19)  and  (20)  into 
(18)  we  find  after  carrying  out  the  integration  over  only  the  term 

with  m = 0 in  sq.  (19)  contributes,  and  obtain 


^ - P^V^N(Ej.)  , (3  co.^  » - 1) 


UD 


where 


00 

= J 4K  3 cos^@-  1)  r'\dV  = (C  - A)  ^ r g(r)  dr  (22) 


is  the  quadrupole  moment  of  the  electrons  near  the  Fermi  level.  The  inter- 
action with  a nuclear  quadrupole  moment  would  involve  a different  quantity  q\ 
in  which  the  charges  of  all  electrons  and  nuclei  are  involved  and  not  just  the 
charges  of  the  conduction  electrons  near  the  Fermi  level. 

The  shift  eq.  (21),  which  averages  to  zero  for  a polycrystalline 
material,  is,  of  course,  superimposed  on  the  isotropic  shift  given 

by  eq,  (3). 

It  is  to  be  noted  that  the  angular  dependence  of  eq.  (21)  is  the  same, 
as  coiUd  be  inferred  in  a purely  phenomenological  way,  if  one  assumes  that 
a linear  tensor  relation  exists  between  the  internal  field  hH  and  the  external 
field  H^.  Since  in  first-order  approximation  only  the  component  of 
parallel  to  will  alter  the  resonance  frequency  we  can  write  down  immediately 

2 . 2 
^ ^ cos  0 ^ *in  0 


(23) 
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The  phenomenological  quantities  can  be  identified  as 

^11  “ 1^-  + 2B^V  N(E_.  q » 

11  18  ^ o ' F)  ^ o 


IS 


o 


9 


The  number  of  nuclei  in  the  powdered  specimen,  for  which  the  magnetic  field 
makes  an  angle  Q ~ith  the  crystallographic  i-ajcis  i»,  of  course,  proportional 
to  sin  0 . This  distribution  over  angles  must  be  transformed  into  a distribu- 
tion function  over  resonance  frequencies.  The  straight-forward  procedure  is 

27 

the  same  as  used,  e.g.,  by  Pake  for  the  shape  of  the  resonance  of  the  two- 
proton  systems  in  powders.  Using  eq.  (23)  we  find  for  the  line  intensity 


2/ (v)cc  sinO 


1/2 


it 


or 

-1/2 

f/(w)  ac(Asj  + - v)  if  j > » - j j 

^ (y)  =0  outside  this  interval. 


(24) 


The  line  shape  Vi*)  is  given  by  the  dotted  line  in  Figure  7 for  the  case 

corresponding  to  a positive  value  of  q.  Dipolar  broadening  and 
finite  life'^time  will  change  the  distribution  to  one  indicated  by  the  solid 
line.  This  is  in  excellent  agreement  with  the  line  shape  observed  in  tetra- 
gonal tin.  We  conclude  therefore  that  q is  positive.  The  dominant  contribu- 
tion to  the  gradient  of  the  electric  field  at  the  nucleus  comes  from  electrons 

28 

along  the  tetragonal  axis.  This  agrees  with  the  crystal  structure  where  the 
four  nearest  neighbors  form  a tetrahedron,  squashed  in  the  direction  of  the 
tetragonal  axis,  while  the  two  next  nearest  neighbors  are  in  the  direction  of 
the  tetragonal  axis.  If  one  assumes  a uniform  charge  distribution  around 
each  nucleus,  it  is  elongated  along  the  tetragonal  axis.  The  experiment  gives 
for  the  anisotropy 


11 


o 


- 0.05% 
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21  29 

We  calculate  N(E_)  from  the  electronic  specific  heat  data  ’ with  tne 
* 12-1 

result  that  V N(E„)  = 0.  19  x 10  erg  Substituting  these  values  into 
° - 24-3 

eq.  (23),  we  obtain  q = + 8.8  x 10  cm  The  relative  anisotropy  of  the 
shift  amounts  to  7 per  cent.  This  should  be  considered 

as  large.  The  p-character  of  the  wave  function  must  be  appreciable.  The 
hyperfine  interaction  in  the  p-state  is  smaller  than  in  the  s-state,  and  further- 
more we  detect  only  the  anisotropy  in  the  p-wave  functions.  The  total  amount 
bf  p-character  should  well  exceed  50  per  cent  near  the  Fermi  level  in  white 
tin. 

In  hexagonal  structures  where  the  nearest  neighbors  have  a nearly 

cubic  arrangement,  much  smaller  anisotropies  should  be  expected,  and  in 

fact  no  other  asymmetric  Lines  have  been  reported.  The  asymmetry  should 

become  more  pronounced  in  higher  fields.  The  problem  treated  here  is 

somewhat  similar  to  the  dipolar  interaction  in  solid  hydrogen  treated  by 
, 30 

Reif.  For  the  proton-spin  pair  a -symmetric  line  always  results  because 
the  spin  of  the  proton  has  equal  probability  for  parallel  or  antiparallel 
orientation.  In  our  problem  of  the  nuclear  spin-electron  spin-pair  inter- 
action, we  have  to  take  the  average  value  of  the  electron  spin  moment,  as 
the  motion  of  the  electrons  is  much  more  rapid  than  the  nuclear  {.recession 
frequency.  The  net  moment  is  always  oriented  parallel  to  This  leads 

to  an  asymmetric  broadening,  proportional  to  the  value  of  H^.  These  two 
features  enable  us  to  distinguish  the  anisotropic  Knight  shift  from  quad- 
rupole  interactions,  in  case  the  nuclear  spin  If  i*  interesting  that 

a value  q for  the  upper  conduction  electrons  can  be  derived  from  this  effect, 
but  this  is  of  no  direct  avail  in  the  absolute  determination  of  nuclear  quad- 
rupole  moments.  As  pointed  out  before,  another  q involving  the  total  charge 
distribution  is  effective  in  that  case. 

The  derivation  for  the  anisotropy  of  the  Knight  shift  showed  no 
relation  with  the  anisotropy  in  the  observed  magnetic  susceptibility.  The 
latter  is  mostly  due  to  diamagnetic  effects  of  the  orbital  motion.  It  is 
possible  that  some  anisotropy  in  the  Pauli  spin  paramagnetism  occurs 
because  of  an  anisotropic  g-factor.  If  the  spin-orbit  coupling  dt  e to  the 
partial  "non-s  character"  in  the  wave  function  is  appreciable  compared  to 
crystalline  field  splittings  in  the  conduction  band,  one  will  have  a g-factor 
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different  from  2,  and  in  a tetragonal  crystal  one  has  to  distinguish  between 
g^j  and  gp  The  partial  s character  of  the  wave  function  then  leads, 
accordin^to  eq.  (3),  to  an  anisotropy  in  the  Knight  shift 


C A V^N(E) 


(24a) 


This  anisotropy  effect  should  be  added  algebraically  to  that  given  by  eq.  (21). 
Phenomenologically  the  two  effects  are  indistinguishable.  It  'viH  be  very 
difficult  to  obtain  experimental  information  about  g^^  and  g^.  In  the  absence 
of  a theoretical  evaluation  we  have  used  eq.  (21)  only  to  obtain  a value  for 
q.  This  will  be  incorrect,  if  the  anisotropy  in  g is  several  per  cent.  This 
would  still  require  a sizeable  amount  of  "non-s  character"  of  the  wave 
function,  but  perhaps  less  than  the  use  of  eq.  (21)  would  indicate.  The 
anisotropic  diamagnetic  susceptibility  will  of  course  give  rise  to  an  aniso- 
tropic magnetic  shielding  of  the  nuclei.  But  the  order  of  magnitude  of  this 

31 

effect  is  much  smaller,  perhaps  a few  parts  in  a million. 

We  finally  extend  the  preceding  considerations  to  the  general  case 
when  there  is  no  axis  of  symmetry.  Instead  of  eq.(20)  we  now  have  for  the 
average  density  of  the  conduction  electrons  near  the  Fermi  level 


=g(r)(Ax^  + By^  + Cz^)  (25) 

Substituting  eqs.  (19)  and  (25)  into  eq.  (18)  one  finds  that  is  now 

not  only  a function  ofO,  but  also  of  4. 


2 A • -jB)  j^TT  (3  cos ^0-1)  +^ir(A-B)sin^0  cos  2^ 


I r g( 


r)  dr 


This  expression  reduces  properly  to  eq.  (21)  for  A ~ B. 

In  order  to  calculate  the  line  shape  we  introduce  three  frequencies 
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Vp  V2  and  v^,  at  which  the  resonance  occurs,  if  points  along  one  of 
the  axes.  We  assume  without  loss  of  generality  that  The 

resonance  frequency  for  the  direction  0,  i is  then  given. by 

V = V j sin^  0 sin^  d + V2  ® ^ ® (27) 

A comparison  with  eq.  (26)  readily  relates  these  phenomenological  quantities 
with  the  A,  B and  C and  with  the  isotropic  shift  • We  first  calculate 
from  eq.  (27)  what  the  Une  shape  would  be  if  0 were  held  fixed 


for 


-1/2  -1/2 

2 2 2 2 
>^cos  0-VjSin  0)  (-v+v^cos  O + v^^in  0) 


2 

> cos 


if  vciv  j 


and 


>cos 


2 


0>0 


if  V v<V2 


and  zero  outside  these  intervals  for  0.  To  find  the  actual  shape  we  have  to 
integrate  over  the  permissible  range  of  0 values 


^ ~ J O' 


0 dO 


and  find  after  some  manipulation  that  !T(  v)  can  be  brought  into  the  form  of 
the  complete  elliptic  integrals. 


U V - ^ ’(Vj-V^)  (v  - Vj)  ^ 


1/2 


dy  for 

(28) 


(v)«e 


f(v3-vp(v3-v2) 
r^3~v  ) vj) 


IT  / 2, 


( (v-Vj)(v3~V2)  2 

^V3-v)  (v^  - Vj)  ^ 


1/2 


dy  for  Vj<  v< 
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^ (v)  = 0 for  vCVj 

v>»3 

The  line  shape  ha*  been  plotted  in  Figure  8,  for  the  case  3 " 2‘ 

3 

The  values  of  the  elliptic  integrals  can  be  found  in  the  tables  of  Jahnke  — Emde. 
There  is  a logarithmic  divergence  at  the  frequency  Further  broadening 

of  the  line  by  dipolar  interaction  will  smooth  this  out  almost  completely, 
and  the  real  line  shape  is  indicated  by  the  drawn  line.  Note  that  the  line  would 
be  symmetrical,  if  ^3  ~ ^2  ~ ^2  ~ ^1*  observations  on  a metal  with  less  than 
ajcial  symmetry  have  been  made.  It  would  be  interesting  if  such  a line  shape 
could  be  identified. 

5 . The  Nuclear  Resonance  in  Thallium  Alloys 

After  the  line  in  metallic  thallium  had  been  found,  experiments  on 
thallium  alloys  were  immediately  successful.  The  lines  retained  their  intensity, 
as  quadrupole  effects  are  absent.  The  width  varied  somewhat  in  the  different 
alloys,  but  not  by  spectacular  amounts.  An  exception  is  perhaps  formed  by 
the  thallium-bismuth  alloys  where  the  line  width  increased  by  a factor  3 to 
110  kc/s.  A remarkable  result  was  obtained  in  a liquid  alloy  (90%  Hg,  10%  Tt) 
at  r.  om  temperature.  The  width  has  the  same  order  of  magnitude  as  in  the 
solid  phases  and  is  apparently  not  narrowed  by  the  rapid  motion  in  the  liquid. 
This  is  all  the  more  remarkable  a*  a spin  lattice  relaxation  mechanism  is 
again  ruled  out  by  the  lack  of  variation  with  temperature. 

Data  on  the  width  and  Knight  shift  are  assembled  in  Table  11.  The 
Knight  shift  exhibits  appreciable  variations  on  alloying.  In  one  phase  it 
changes  continuously  with  composition.  This  behavior  is  illustrated  for  the 
thallium-tin  system  in  Figure  9.  The  solubility  of  thallium  in  tin  is  too 
small  to  obtain  data  at  the  other  end  of  the  phase  diagram.  The  relative 
shift  increases  from  1.54  per  cent  in  pure  thallium  to  1.97  per  cent  at 
the  phase  boundary.  The  Knight  shift  of  tin  in  the  alloy  (20%  Sn-80%  Tt ) on 
the  other  hand  is  smaller  (0.67%)  than  in  pure  tin  (0.75%).  According  to 
eq.  (3)  the  variation  can  be  caused  by  a change  in  the  quantities  C,  and 
N(E).  In  a first  approximation  it  is  probably  correct  to  assume  that  CV^ 
is  constant,  as  the  hyperfine  interaction  normalized  per  atom  is  inversely 
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proportional  to  the  atomic  volume  V^.  Then  the  plot  of  the  Knight  shift 

is  also  a plot  of  density  of  states  near  the  Fermi  level  N(E_).  For  a free 

1/2  * 

electron  gas  this  quantity  is  proportional  to  E_  . It  is  plausible  that  the 

X 

addition  of  tin  to  thallium  raises  the  number  of  conduction  electrons  per 
unit  volume  and  consequently  Ep  andN(Ep).  There  are  some  data  onlhe  mag- 
netic susceptibility  which  also  indicate  an  increase  of  N(E).  The  specific  sus- 

» =-  0.  llx  10'^. 
Sn 

The  density  N(E)  levels  off  near  the  phase  boundary;  a decrease  which  might 
be  expected  according  to  Jones'  explanation  of  the  Hume-Rothery  rules  is 
not  observed.  Additions  of  less  than  3 per  cent  tin  do  not  change  the  Knight 
shift.  The  density  of  states  N(Ep)  remains  constant  at  first  and  increases 

only  when  more  than  3 per  cent  tin  is  added.  This  is  in  agreement*  witli 

34 

1;  riedel's  theory  of  impurities  in  metals.  The  tin  atoms  in  small  concentra- 
tion act  like  impurities,  which  have  either  a bound  electron  or  a compressed 
Fermi  gas  of  conduction  electrons  around  them.  In  either  case  the  maximum 
momentum  and  therefore  N(Ep)^are  unchanged.  For  higher  concentrations 
the  impurities  interact  and  N(Ep)  increases. 

The  line  width  increases  by  a factor  two  if  10  per  cent  tin  is  added.  Data 
are  given  in  Figure  10.  The  thallium  nuclei  in  the  solid  solution  do  not  all 
have  the  same  environment.  There  will  be  a certain  distribution  in  the  value 
of  |i|<<0)|^  depending  on  the  presence  of  tin  atoms  in  the  immediate  neighborhood. 
This  will  lead  to  a distribution  of  Knight  shifts.  Although  one  would  expect 
an  asymmetric  distribution  for  small  tin  concentrations,  an  asymmetric  broad- 
ening has  not  been  observed. 

Discontinuities  in  the  Knight  shift  occur  when  phase  transitions  take 
place.  In  the  two -phase  regions  of  the  phase  diagram  the  resonance  will 
therefore  be  double.  This  is  illustrated  in  Figure  II,  where  the  Knight 
shift  is  plotted  versus  composition  of  thallium-mercury  alloys.  The  rela- 
tive change  in  the  Knight  shift  is  more  than  70  per  cent  in  going  from  the  p - 
to  the  y - phase.  The  change  is  partly  caused  by  volume  effects  but  mostly 
by  electronic  structure.  Another  example  of  the  effect  of  phase  change  on  the 
Knij;ht  shift  is  the  melting  of  T6  - Hg  alloys.  A pure  increase  in  volume  would 

Professor  Brooks  first  proposed  this  explanation. 


ceptibility  increases  from  = - 0. 22x  10  ^°^g0%Tt  -t  20% 
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Table  U 


Phase 

or 


System 

Composition 

Isotope 

Knight  Shift 

Na-Te 

pure  T6 

T(205 

1.54^.05 

NaTt 

„205 

i 

-1.03 

Na^Te 

Te^OS 

.68 

! 

NaTfc 

N.“ 

- .016 

iMg-Tt 

MgTt 

2.  12 

Iln-Te 

50  At.%  In-Te 

1.68 

;Sn-Tt 

0-25  At.%  Sn-Te 

1.54  to  1.97 

see  det 

ailed  shift  vs . 

composition  data 

25  At.%  Sn-Te 

s„“’ 

. 67 

pure  p-Sn 

s„“’ 

.74 

pure  a Sn 

s„“’ 

0.0 

Hg-Te 

13  At.%  Hg-Tt 

1.71 

60  At.%  Hg-Te 

Te  Viquid) 

1.49 

72  At.%  Hg-Te 

1.01 

92  At.%  Hg-Te 

T£  ^^^(liquid) 

1.18 

Pb-Tt 

34  At.%  Pb-Te 

Tj205 

1.39 

90  At.%  Pb-Te 

Te“* 

1.90 

pure  Pb 

Pb“’ 

1.24 

Bi-Te 

6 At.%  Bi-Tt 

TjZOS 

1.71 

19  At.%  Bi-Te 

1,74 

59  At.%  Bi-Te 

Tt“5 

1. 1 

between  points  of  maximum  slope 


second  moment 


PERCENT 


PHASE  COUNDARY> 


ATOMIC  PERCENT  THALLIUM  IN  TIN 


Figure  9 The  shift  oC  the  resonance  in  solid  solutions  of  tin  in 

thallium.  The  fact  that  tin  concentrations  less  than  3%  have  no  effect 
on  the  Knight  shift  can  be  considered  as  proof  of  Friedel's  theory  of  im- 
purity levels  in  metals,  which  leave  the  Fermi  level  unaffected.  For 
higher  concentrations  the  curve  represents  the  filling  up  of  the  conduction 
band  and  gives  the  density  of  electron  states  N(E). 
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presumably  lower  the  Knight  shift  according  to  eq.  (3).  Actually,  an  increase 
of  20  per  cent  is  observed  on  melting,  indicating  a change  in  the  electronic 
structure.  The  behavior  of  the  Knight  shift  over  the  phase  diagram  is  roughly 
what  could  be  expected.  The  knowledge  of  the  electronic  structure  of  these 
alloys  is,  however,  insufficient  to  attempt  a quantitative  explanation. 

A thallium-magnesium  alloy  c>.hibits  the  largest  Knight  shift  of  2.1  per  cent. 

A most  remarkable  result  was  obtained  in  the  ordered  Na  Tt  structure. 

This  intermetallic  compound  shows  a negative  Knight  shift  ( - 1.03%)  for  thallium. 

23 

Since  it  is  an  ordered  structure  with  cubic  symmetry,  the  Na  resonance 
could  also  be  observed,  and  it  shows  a small  shift  towards  lower  frequencies. 

A shift  of  - 0.016  per  cent  was  found  compared  to  a solution  of  Na  I.  In 
metallic  sodium  an  ordinary  positive  shift  of  0.  1 1 per  cent  occurs.  We  hav» 
not  yet  performed  experiments  at  other  frequencies  and  do  not  know  whether 
the  shift  is  really  proportional  to  v^.  although  it  seems  unlikely  that  a chemical 
shift  would  be  this  large.  Experiir  :nts  on  other  alloys  with  the  Na  Tt  structure 
are  also  planned. 

A possible  explanation  of  the  negative  shift  might  perhaps  be  found  along 

35  -f-f 

the  following  lines.  Abragam  and  Pryce  have  shown  that  Mn  ions,  and 
other  ions  of  the  iron  group,  exhibit  configurational  interaction,  in  which  a 
3s-eiectron  is  promoted  to  the  4-s  orbital.  The  total  spin  or  the  magnetic 
moment  of  the  ion  is  not  changed  by  this  promotion,  but  the  interaction  with 
the  nucleus  is  changed.  The  electron  spin  interacts  more  strongly  in  the 
3 8 orbit  than  in  the  4-s  orbit.  The  effective  spin  density  at  the  nucleus 
due  to  the  unpaired  s electrons  is  determined  by  overlap  integrals  of  radial 
wave  functions.  This  configurational  interaction  accounts  for  the  hyper- 
fine  structure  of  the  Mn^^  -ion.  Abragam  and  Pryce  show  that  it  is 
proportional  to  the  amplitude  of  the  admixed  wave  function  and  not  to  its 
square.  It  lo  an  interference  effect  which  actually  gives  a negative  sign 
for  the  A I • interaction  in  the  Mn  ion.  The  question  is  whether  con- 
figurational interaction  in  a metal  might  cause  some  mixing  with  a state 
in  which  an  inner  s electron  is  promoted  to  higher  orbits  in  the  conduction 
band  with  s- wave  character,  and  whether  this  effect  would  yield  a particularly 

i 

The  authors  are  indebted  to  Dr.  Abragam  for  calling  their  attention  to  this 
effect  and  for  making  the  manuscript  available  before  publication. 
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large  negative  value  for  Na  Tt  . Clearly  much  theoretical  work  and  more 

experimental  data  are  needed  to  confirm  this  hypothesis.  The  width  of 
23 

the  Na  resonance  in  Na  Tt  is  normal.  It  should  be  reiitembered  that  the 
tabulated  experimental  and  theoretical  values  are  not  cirectly  comparable. 
The  width  of  the  thallium  resonance  is,  however,  of  the  same  order  of 
magnitude  as  found  in  all  other  samples.  This  is  another  argument  that  the 
anomalously  large  width  is  a characteristic  of  the  atomic  configuration  of 
the  heavy  elements,  rather  than  of  the  electronic  structure  of  the  conduction 
band. 

Unfortunately,  we  have  no  data  on  the  spin -lattice  relaxation  time  for 
any  ol  the  observed  resonances.  Wo  check  on  the  validity  of  the  Korringa 
relation  is  therefore  possible. 


The  authors  wish  to  express  their  indebtedness  to  Professors  H.  Brooks, 
J.  Friedel,  H.  Jones  and  Dr.  C.  Herring  for  several  valuable  discussions. 
One  of  the  authors  (T.  J.R.)  acknowledges  a stipend  from  the  California 
Research  Corporation. 
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